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ABSTRACT 

The high-magnification gravitational lensing cross-section of a galaxy is enhanced sig- 
nificantly if a baryonic disc is embedded in its dark matter halo. We investigate the 
effects of a population of such discs on the probability of detecting strongly lensed 
images of distant galaxies and quasars. The optical depth to lensing is always more 
than doubled. The effects are particularly significant at magnifications greater than 
about 50, for which the optical depth typically increases by a factor greater than 
5. If either the fraction of discs or the typical disc-to- halo mass ratio increases with 
redshift, then the optical depth is expected to be increased even further. Obscura- 
tion by dust in the lensing disc is expected to counteract the associated enhanced 
magnification bias for surveys in the optical waveband, but not in the radio and mil- 
limetre/submillimetre wavebands. The presence of galactic discs should hence lead to 
a significant increase in the number of lensed galaxies and quasars detected in both 
radio- and millimetre/submillimetre-selected surveys. An increase by a factor of about 
3 is expected, enhancing the strong case for millimctrc/submillimetre-wave lens sur- 
veys using the next generation of ground-based telescopes and the FIRST and Planck 
Surveyor space missions. 

Key words: galaxies: evolution - galaxies: fundamental parameters - cosmology: 
theory - gravitational lensing - radio continuum: galaxies 



1 INTRODUCTION 

There has been much recent interest in the effects of galactic 
discs on the lensing cross-sections of galaxies (Mailer, Flores 
& Primack 1997; Wang & Turner 1997; Bartelmann & Loeb 
1998; Keeton & Kochanek 1998; Koopmans, De Bruyn & 
Jackson 1998; Moller & Blain 1998). A gravitationally dom- 
inant baryonic disc in the core of a dark matter halo, edge- 
on to the line of sight to a distant source, can significantly 
increase the mass surface density, and thus the lensing effi- 
ciency. 

The image geometry produced by an individual lens can 
be modified significantly by the presence of an edge-on or 
almost edge-on disc; however, the cross-section to multiple 
imaging is only increased by a factor of about 50 per cent 
when averaged over all disc inclinations. In contrast, the 
cross-section to magnifications greater than 10 is expected 
to be increased by a factor of a few after the same aver- 
aging process (Moller & Blain 1998). Hence, depending on 
the abundance and form of evolution of discs in dark mat- 
ter haloes, the strong lensing optical depth at moderate and 
large redshifts could be increased very significantly. Bartel- 
mann & Loeb (1998) have also recently discussed this effect. 
They note that it will probably not be possible to exploit 
any enhanced magnification in the optical waveband because 



the flux densities of the lensed images will be attenuated by 
dust associated with the disc. For the same reason, Mailer 
et al. (1997) did not consider the effect of disc lensing on 
the magnifications of multiple images because of possible 
differential extinction between the images. 

Dust extinction is important in the optical waveband; 
however, a dusty lensing disc is optically thin in the radio 
and millimetre/submillimetre wavebands. Indeed, because 
of fc-correction effects (Blain & Longair 1993), a lensing ob- 
ject is often expected to be fainter than the lensed images 
in the submillimetre waveband (Blain 1997a, 1998b). Hence, 
any enhanced magnification bias due to the effects of discs 
would be easy to detect in these wavebands. Existing ra- 
dio lens surveys, such as the CLASS survey (Browne et al. 
1998), appear to show an overabundance of lens candidates 
that can be identified as highly inclined discs. The potential 
importance of submillimetre- wave lens surveys has been dis- 
cussed for lensing by both clusters (Blain 1997a, 1998a) and 
galaxies (Blain 1996, 1997b, 1998b,c). Submillimetre- wave 
lensing by clusters has now been demonstrated to be a very 
useful phenomenon (Smail, Ivison & Blain 1997). The most 
striking feature of submillimetre- wave galaxy lens surveys is 
that the fraction of strongly lensed galaxies detected in a 
carefully selected sample can exceed 5 per cent (Blain 1996, 
1998b). 
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In Section 2 we discuss the statistical features of lens- 
ing by galaxies containing a disc, using techniques and mod- 
els of lensing by disc galaxies developed by Moller (1996), 
Mailer et al. (1997) and Moller & Blain (1998), and the 
statistical lensing formalism (Peacock 1982; Pei 1995) ap- 
plied in the submillimetre waveband by Blain (1996). In 
Section 3 we predict the effects of a potentially evolving pop- 
ulation of discs (Lilly et al. 1998; Mo, Mao & White 1998; 
Mao, Mo & White 1998; Marzke et al. 1998) on the lensing 
cross-sections. In Section 4 we estimate the surface density 
of strongly lensed galaxies and quasars in the light of recent 
submillimetre- wave observations by Smail et al. (1997) us- 
ing the new SCUBA camera at the James Clerk Maxwell 
Telescope (Holland et al. 1999). See Smail et al. (1998, 
1999, in preparation), Ivison et al. (1998), Blain, Ivison & 
Smail (1998a), Blain et al. (1998b, 1999a,b) and Frayer et al. 
(1998) for a discussion of the latest results. In Section 5 the 
prospects for detecting such images in future surveys using 
forthcoming ground-based and space-borne instruments are 
discussed. 

The values of cosmological parameters have a significant 
effect on the properties of galaxy-galaxy submillimetre- wave 
lensing (Blain 1998a). Unless otherwise stated, we assume an 
Einstein-de Sitter world model with Hubble's constant Ho = 
60kms~ 1 Mpc -1 , in accordance with the models presented 
by Mailer et al. (1997) and Moller & Blain (1998). Note 
that the abundance of submillimetre-wave lensed galaxies 
predicted in an Einstein-de Sitter model is less than that 
expected in a model with a smaller density parameter (Blain 
1998a). 



2 THE STRONG LENSING OPTICAL DEPTH 
2.1 Basic formalism 

Millimetre/submillimetre-wave source counts of lensed im- 
ages were predicted by Blain (1996), using formalism de- 
veloped by Peacock (1982) and Pei (1995) to estimate the 
optical depth to strong lensing attributable to galaxies. In 
this formalism the probability that a source at redshift z is 
lensed by a magnification between /i and [i + d[i is F(/i, z). 
This function is obtained by integrating the cross-sections of 
individual lenses to magnifications greater than /i, o>(M, z) 
in the plane of the lens, both along the line of sight to z and 
over the mass function of lenses N, and then differentiating 
with respect to (i, 

dr 



a fi (M,z')(l + z'YdN(M) — dz'. (1) 



d/i 

Here r is the comoving radial distance element. The number 
and flux density conservation conditions, 



/ F(p, z) d/i 
Jo 
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Jo 



nF([i,z) d/i = 1, 



(2) 



must be satisfied. If the intrinsic luminosity function of 
galaxies $(1/, z) is known, then the source count of unlensed 
galaxies 

L m in is the luminosity of a source at a redshift z that pro- 
duces a flux density S v ; zo is the maximum redshift of the 



source population; and D(z a , z D ) is the angular diameter dis- 
tance between redshifts z a and Zb- The count of lensed galax- 
ies can be determined using equation (3) , if $ is replaced by 
an effective luminosity function, 



*'(L,z) = 



d/t, 



(4) 



corrected to include the effects of lensing. p m ax, the maxi- 
mum magnification, is determined by the intrinsic size of the 
source and the geometry of the lens configuration. In gen- 
eral, larger values of /i max are expected for smaller sources. 
A*max(z = 1) — 25 is expected for lensing by a singular 
isothermal sphere (SIS) lens if the intrinsic size of the source 
is lkpc (Peacock 1982). 



2.2 Lensing by singular isothermal spheres (SISs) 

The cross-section to high-magnification lensing by an SIS 
at redshift z' for a source at redshift z, as measured in the 
plane of the lens, 



af S (/i,M ) =64^ 
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(5) 



(Schneider, Ehlers & Falco 1992). M e is the mass of the lens 
enclosed within a radius R. The distance ratio inside the sec- 
ond square bracket is referred to as Dr(z', z) in subsequent 
equations. In the case of lensing by SISs, the magnification 
and redshift dependences of F (equation 1) can be separated 
at large magnifications, yielding the strong lensing probabil- 
ity -Fsis = a(z)/i -3 that is applicable here. The form of Fsis 
at all magnifications can be approximated by 



if nt < ii, < n 11 
otherwise. 



(6) 



The first term describes the typically small magnification 
correction that must be applied to most sources in order to 
satisfy the number and flux density conservation conditions 
(equation 2). 

The form of the function a(z) depends on both the val- 
ues of cosmological parameters and the form of evolution of 
lensing galaxies (Blain 1996, 1998a). If the distribution of 
lenses in mass M is described by a Press-Schechter function 
(Press & Schechter 1974), 



dN(M, z) = 



h 7 



/ M \ 7/2 r ( M\ 

(w) exp [- 



dM, (7) 



then equation (1) can be evaluated analytically. M*(z) in- 
cludes all the details of the evolution of structure resulting 
from the growth of perturbations in an expanding Universe; 
7 = 4/3 if the primordial density fluctuations are scale- 
independent; p g is the mean smoothed density of bound ob- 
jects in the Universe; and M is the mass of a region of the 
Universe that has turned around from the Hubble flow and 
is in the process of collapse at a particular epoch. M is re- 
lated to the mass enclosed within a radius 7?, M e (equation 
5) by the constant e, that is, 



M = eM. 



(8) 



M* is the equivalently transformed version of M* . This 
transformation allows the masses in equations (5) and (7) to 
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be linked in the same expression. For example, if the pop- 
ulation of lensing galaxies consists entirely of SISs, then by 
evaluating equation (1), 



a(z) = 128\/7rpge 
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Dl{z' \z)(l + z 'f M :(z')^dz' . 



(9) 



If -Fsis is given by equation (6) then 

*s IS = - $ ( -) + «W f maXW ^ * (-) d/x, (10) 

in which /(//) = /i~ 3 . This equation was used to derive the 
results for the probabilities of statistical lensing presented 
by Blain (1996, 1998a). 



2.3 Lensing by discs 

Moller & Blain (1998) recently presented the results of simu- 
lations of the high-magnification lensing cross-section for an 
exponential disc embedded within an SIS halo. When the 
disc was within about 30° of being edge-on to the line of 
sight, this cross-section was found to be increased by a sig- 
nificant factor as compared with that for an SIS halo alone. 
An example of the relative magnification a^/a^ 5, is plotted 
in fig. 3(a) of Moller & Blain (1998), for a disc/halo combi- 
nation well matched to that of the Milky Way. We have now 
extended these simulations to include a wide range of halo 
and disc masses, and so to investigate the statistical prop- 
erties of lensing by such objects. The central surface mass 
density of the disc within the halo was kept constant (Mao 
ct al. 1998), while the scalelength of the disc was varied in 
such a way that the ratio of the mass of dark matter in the 
SIS halo within a radius R = lOkpc and the mass of the 
disc remained constant. The results are now parametrized 
to generate an analytical description of the effects of lensing 
by a population of discs. 



2. 3. 1 Formalism 

The relative cross-section to lensing by a magnification 
greater than \i due to a combination of a disc and an SIS 
halo, as compared with an SIS halo alone, was approximated 
by a function of the form 



{(3 - a)n. 



1+ V 



3+a 



-£ S = C(M C )< 
where 

C(M e ) = K 



1 + a + (3 + a/3 1 + a ' 



a- 13 3+/3 



1 + 



Me 

Mo 



if fj, < 
otherwise, 



(11) 



(12) 



Any systematic evolution of the scalelength/mass of discs 
within dark matter haloes can be incorporated by modifying 
the value of Mo, defined as the mass of dark matter enclosed 
within lOkpc of the core of an SIS halo, as a function of 
redshift. 

Equation (11) has the useful property that when dif- 
ferentiated to obtain da^/djj, (required in equation 1), the 
result has a simple double-power-law form, 



Table 1. The values of the best-fitting parameters derived from 
the extended results of Moller & Blain (1998). The associated 
best-fitting curves are shown in Fig. 1. The best-fitting values of 
the worst constrained parameters, (3 and fit,, are related by the 
equation /3 = 0.017/i* - 2.6. 



Parameter 


Value 


Parameter 


Value 


a 
V 

M 

M* 


-2.85 ±0.16 
-0.87 ±0.06 
(2.6 ±0.4) x 10 11 Mq 
2.9 x 10 11 M 


P 

A** 

K 


9 19 +0.25 

1.7 ±0.5 



do> 
d/i 



64tt 



G M c 
1?~R 



' C(M C ) | 



if n < n* ; 
otherwise. 



(13) 



The parameters a, (3, ju*, K, Mo and r\ were determined 
by fitting the function in equation (11) to the results of the 
numerical simulations (Moller & Blain 1998). The form of 
the best-fitting models, and probability contours determined 
by fitting to the worst-constrained pairs of parameters, (3 
and /it*, and r\ and Mo, are all shown in Fig. 1. The values 
of all the best-fitting parameters are listed in Table 1. 

The resulting model for the cross-section to strong lens- 
ing by discs within SIS haloes can be used to estimate the 
optical depth to strong lensing Fdi sc . Fdisc is obtained as a 
function of redshift by evaluating equation (1), incorporat- 
ing the fraction of haloes that contain a baryonic disc f&{z): 



Fdisc — — 



d 



d/i _ 

[(l-/d)< S (M) + / d <r M (M)] dN(M)dz'. (14) 
For magnifications /i > fi t this can be written as, 

Fdisc = [ai(z) - a 2 (z)] I n (fi) + [a 2 {z) + a 3 (z)] Fi(m)> ( i5 ) 

the sum of two magnification-dependent terms, I n and Id, 
which are modulated by redshift-dependent terms that are 
combinations of three functions, ai, 0,2 and 0,3. As in equa- 
tion (6), a low-magnification correction term containing H 
and no must be included in order to approximate the com- 
plete probability function. 

I n and Id are associated with lensing by pure SIS haloes 
and by discs within SIS haloes respectively: 
-3 



UP) 

and 

w=f {£-v 

The a-functions (equation 15) are 



if (j, < fi* ; 
otherwise. 



ai(z) =r 

a 2 (z) = r 
and 

a 3 (z) = r 



1 1 

7 + 2 

1 1 

7 + 2 



77 + I 



J{z')M*{z')dz, 



f J(z')fd(z')M:(z')dz', 
Jo 



7 2 
in all three of which 

J(z') = 128v^Fpe 



M*(z') (1+v) 

, / x fd(z')dz', 

y ' M {z')v 1 v 1 



' Dl{z>,z){l + z'?^; 



(16) 
(17) 

(18) 
(19) 

(20) 
(21) 
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Figure 1. The best-fitting models that describe (a) the magnification dependence of the function o"^i/c^ IS (equation 11) and (b) 
the mass dependence of the function C(M e ) (equation 12), shown by the smooth curves. The results of the simulations arc shown 
by the jagged line in (a) and the empty squares in (b) respectively. The quality of the constraints imposed on the worst constrained 
pairs of parameters in each case are shown in (c) and (d) respectively. The three contour levels correspond to likelihoods la, 2a 
and 3(7 away from the maximum value, which is marked by a cross. The values of the best-fitting parameters are listed in Table 1. 



Note that if discs do not evolve, then the redshift depen- 
dences of ai , a 2 and a 3 are identical. 

The relative normalization of 03 as compared with ai 
and 02 depends on the ratio of the mass Mo, derived from 
the results of the numerical simulations described above, 
and the typical mass in the Press-Schechter function M* 
(equations 7 and 8). The value of Mg can be determined 
by using the Tully-Fisher relation (Hudson et al. 1998), an 
empirical power-law relation between the B-band luminosity 
and circular velocity of galaxies. A circular velocity v c ~ 
360 kms -1 is indicated by the Tully-Fisher relation for an L* 
galaxy with a B-band magnitude Mb = —21 at the present 
epoch. If L* is assumed to be equivalent to M* , then Ml 
can be determined by comparison with the Milky Way, for 
which M C (R < lOkpc) = 10 11 M and v c ~ 210 kms -1 . 
M* ~ 2.9 x 10 11 M Q is expected in this case. 

By combining equations (15) to (21), we can model the 
optical depth to strong lensing by a population of SIS dark 
matter haloes that contain an evolving population of bary- 
onic discs. A new effective luminosity function & disc can be 
produced by including the general form of fdisc in equa- 
tion (4). In order to evaluate the count of lensed sources 



using this function, the forms of 01(2), 02(2), 03(2), fd(z) 
and Mo (2) are required. The redshift dependences of the a- 
functions in both hierarchical and non-hierarchical models, 
and three different world models, are illustrated in Fig. 2. 
The form of evolution of discs, discussed by Lilly et al. (1998) 
and Mao et al. (1998), is considered in Section 3.1. 

We parametrize the fraction of haloes that contain discs 
fd = /o(l + z) p and the scalelength r s oc (1 + z) q/2 . The 
comoving abundance of discs and the mass of the disc within 
a typical halo are thus expected to evolve as M*(2) _1 /d(2) 
and M*(2)(l + 2) 9 respectively, if the central surface density 
of the disc is assumed to be independent of redshift. In an 
Einstein-de Sitter model, M* oc (1+2) -1 would be expected 
because of the hierarchical evolution of bound objects. 



2.3.2 Normalization of the terms in F disc 

The relative size of the two terms in equation (15) depends 
on the parameters fa , p, q and pe. To clarify the dependence 
of the results on the form of evolution, we introduce the 
parameters r and R T , the total probability of lensing by a 
magnification fi > fxt, and the ratio of the probabilities of 
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[)„=!. 0; tl A =0.0 and p = q = Q in all i 



(b) 

Thick: Q =1.0; tl A =0.0 
Medium: n„=0.3; Q A =0.0 
Thin: Q„=0.3: Q,=0.7 



Solid: a 1 (z) 
Dashed: a 3 (z) 
Dotted: a ? (z) 



jV/ evolution and p = q = in all cases 



Redshift 2 



Figure 2. A comparison of the three redshift-dependent a-functions that contribute to F^ lac (equation 15). In (a) the effects of 
including evolution of M* in the Press-Schechter function are shown; a\ and a.2 are affected in the same way, but the effect on a% 
is much less significant. In (b) M* is assumed to evolve, and the effects of three different world models are compared: see also Blain 
(1998a). The functions are normalized to match the low-redshift value of a(z) used in Blain (1996), with a value of fo = 0.15: see 
Section 3.1. 



lensing by a disc as compared with lensing by a halo alone 
respectively: 

{' 00 P CO 

t = [01 - 02] / In{fJ,)dfi+ [02 +03] / Id(tJ-)dfi, (22) 

J fit J Mt 

and 



"2 



+ a 3 JT Jd (^) d ^ 



ai — a2 J°° J n (At) d/i 



(23) 



At a low redshift z\ , for which all three a-functions are pro- 
portional to z 3 , and for values of fit — 2 then 



G 



c 2 R 



(l + /oT)T 
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Df\( 2 ', 2 )(l + 2 ') 2 M*(z')^7d Z ', (24) 



in which 



T = 



and 



AC 
Mo 



/o 



l-/o 



(1+T). 



(25) 



(26) 



Based on the values of parameters listed in Table 1, T is 
expected to take the value 1.6. If r ~ 10~ J at z = 0.1 
(Blain 1996), then the density parameter associated with p s 



is fLe ~ 8x 10 . Hence, since fL 



8xl0 _/ . This 



would indicate that 8 per cent of the mass of the dark matter 
in the SIS halo of a M* galaxy at the present epoch was 
enclosed within 10 kpc of the core, equivalent to a truncation 
radius of 125 kpc for the SIS. 

If fo = 0.6 (see Section 3.1), then R T ~ 3.6, and so 
about 80 per cent of lensing events should be caused by disc 
galaxies, neglecting the effects of both magnification bias 
and evolution of the population of lenses. 

By combining equations (24) to (26), the low-redshift 
normalization of the a-functions can be determined by com- 
parison with the normalization of a(z) in a disc-free model 



(equation 6): see Blain (1996) for a summary of the esti- 
mated form of a(z). The normalization of ai is given by 

ai(^i) = ir n (^i) M t- 2 (l + /oT), (27) 

where r n is the probability of lensing in a model without 
discs, to which the predictions in Blain (1996, 1998a, 1998c) 
were normalized. The other two a-functions are normalized 
to 



02(21) = /oai(zi), 
and 

fl 3(zi) = foTai(zi). 



(28) 



(29) 



3 LENSING CROSS-SECTION PREDICTIONS 

3.1 The population and evolution of discs 

Lilly et al. (1998) discussed the morphological properties 
of galaxies detected in the Canada-France Redshift Sur- 
vey (CFRS). To avoid surface brightness selection effects, 
they classified galaxies based on their appearance within 
4ft -1 kpc of the core. They found no evidence for evolution 
in the abundance of spiral galaxies out to redshifts greater 
than 0.5. About 18 per cent of the galaxies with identifi- 
able morphologies in the small-radius sample contained a 
disc. Based on a B-band selected local galaxy survey, unaf- 
fected by surface brightness selection effects, Marzke et al. 
(1998) find that about 64 per cent of local galaxies have a 
disc morphology. Brinchmann (private communication) find 
that about 50 per cent of galaxies selected in the I band are 
classified as containing a disc, either by eye or by machine. 
The fraction does not appear to depend strongly on either 
the redshift or the intrinsic luminosity of the galaxies. The 
zero-redshift abundance of discs in haloes, fo, thus appears 
to be about 0.6. 

Based on the results of theoretical and numerical work 
on the evolution of discs (Mo et al. 1998), Mao et al. (1998) 
argued that the results presented by Lilly et al. (1998) are 
consistent with a comoving abundance of discs that varies as 
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0.1 0.3 0.5 1 3 5 0.1 0.3 0.5 1 3 5 



Redshift 2 Redshift z 

Figure 3. A comparison of the effects of disc evolution on the two redshift-dependent functions that contribute to Fdi sc (equation 
15), ai — a2 and 02 + 03. In (a) five different forms of the abundance of discs /d(z) = /o(l + z) p are compared, with p = — 1, 0, 1, 2 
and 3; the disc-to-halo mass ratio does not evolve. In each case, the value of fo is chosen to satisfy the constraint that the fraction 
of discs found in the CFRS i?d ~ 0.18 (equation 30). In (b) the effects of five different forms of the evolution of Mo oc (1 + z) q are 
included, with q = —2, —1, 0, 1 and 2; a value of p = 1 is assumed in all cases, corresponding to a value of fo = 0.089 if i?d = 0.18. 
The discontinuities in the gradients of the solid curves at high rcdshifts in (a) arc caused by /<j reaching its upper limit. Note that 
the strongly evolving models have values of fo that are much smaller the observed zero-redshift value (Marzke et al. 1998). 



(1 + z) 3 , and a typical scalelength r B oc (l + z) -1 . If M* does 
not evolve, then values of p = and q — correspond to no 
evolution of either the abundance or the typical scalelength 
of discs (Lilly et al. 1998); if M* does evolve, then values 
of p = —1 and q = 1 would be required. If the abundance 
and scalelength of discs evolves as indicated by Mao et al. 
(1998), then values of p = 2 and q = — 1 are required if M* 
evolves; values of p — 3 and q — — 2 are required if M* does 
not evolve. 

Several physical processes can be invoked to explain the 
evolution of discs. If discs grow by the accretion of baryons 
from their haloes, then the disc-to-halo mass ratio would be 
expected to be smaller at large redshifts. Disk-disc or disc- 
halo mergers would generally disrupt discs, reducing the 
abundance of discs at low redshifts. Low-mass discs could 
also be disrupted by supernova-driven winds associated with 
episodes of powerful star formation. A determination of the 
evolution of the properties of discs would thus have impor- 
tant consequences for understanding the formation and evo- 
lution of galaxies. 

The dominant contribution to the lensing cross-section 
is provided by galaxies at redshifts less than unity. Hence, 
the results of Lilly et al. (1998), which sample this redshift 
range, should provide useful information about the normal- 
ization constant fo- In the absence of selection effects, the 
fraction of discs, 

I /l(1+2) (iW M {z) ^ dz 

Rd ^ j , (30) 

subject to the constraint that /o(l + z') p < 1. If p = 0, 
then R d (~ 0.18) = / /(l - fo), and so f ~ 0.15. If the disc 
fraction increases with increasing redshift, that is p > 0, then 
the value of fo will be reduced. The form of evolution of M* 
has very little effect on the results. For values of p = — 1, 0, 
1, 2 and 3, fo takes the values of 0.25, 0.15, 0.089, 0.051 and 
0.025 respectively. 



Table 2. The values of the parameters in the disc lensing models 
used to derive the results in Figs 4 and 5. Both models are con- 
strained by requiring that they provide a reasonable, if schematic, 
account of the abundance of discs determined in the CFRS fields 
by Lilly et al. (1998), normalized to the disc fraction at zero red- 
shift found by Marzke (1998). In both models, discs have the same 
properties at all rcdshifts less than unity (Lilly et al. 1998). See 
Section 2.3.1 for a definition of the parameters fo, p and q. 
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In the absence of universal agreement on the interpreta- 
tion of the CFRS data, we present forms of the functions oi, 
02 and a-i derived for a range of different values of p and q to 
illustrate the possible effects of different forms of disc evolu- 
tion: see Fig. 3. The most significant effects are produced by 
large values of q. These values correspond to models in which 
the disc-to-halo mass ratio evolves strongly. Large lens sur- 
veys carried out in the radio and millimetre/submillimetre 
wavebands should eventually provide a large sample of mass- 
selected lensed objects. These samples could be used to re- 
solve any remaining uncertainty in the interpretation of ex- 
isting surveys. 



3.2 Forms of F d i sc derived for evolving discs 

The optical depths to strong lensing i^isc and Fsis can now 
be calculated and compared. In Fig. 4 Fkisc/^sis is shown 
for both models of disc evolution listed in Table 2 as a func- 
tion of both magnification and redshift. The optical depth 
to strong lensing is increased by incorporating the effect of 
lensing by discs by a factor between 2 and 3 at redshifts of 
about unity and magnifications of several tens. The effect 
of including discs is more significant at greater magnifica- 
tions. The consequences for the counts of lensed galaxies 
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Figure 4. The ratio of the differential optical depths to high-magnification lensing produced by a population of discs within SIS 
haloes and a population of pure SIS haloes, Fdisc/^SIS as a function of redshift and magnification. Models A and B, which include 
different values of the parameters p and q, are defined in Table 2. The lensing parameters listed in Table 1 arc assumed. The results 
are not affected significantly by changes to the world model. 



and quasars in the submillimetre waveband are discussed 
below. 

4 LENSED GALAXIES IN THE 
SUBMILLIMETRE WAVEBAND 

4.1 Introduction 

The selection functions in submillimetre-wave surveys are 
very broad and extend out to high redshifts. This is because 
the flux density-redshift relation for distant dusty galax- 
ies is expected to be fiat over the range of redshifts from 
about 0.5 to 10, owing to large negative fc-corrections. These 
arise when dust emission spectra, peaking in the rest-frame 
far-infrared waveband, are redshifted into the submillime- 
tre waveband (Blain & Longair 1993). A submillimetre-wave 
survey thus selects high-redshift galaxies preferentially, and 
because a flat flux density-redshift relation naturally cor- 
responds to a steep source count, magnification biases in 
such a survey are expected to be large (Blain 1996). At that 
time, however, source counts of submillimetre-wave sources 
had not been determined. Now Smail et al. (1997) have de- 
tected a population of distant dusty galaxies, including a 
dust-shrouded starburst/ AGN at a redshift z = 2.803 (Ivi- 
son et al. 1998; Frayer et al. 1998). A fuller description of 
the survey (Smail et al. 1998, 1999, in preparation; Blain 
et al. 1999b) and its consequences (Blain 1998c; Blain et 
al. 1998a, b, 1999a) can be found elsewhere. The results of 
other submillimetre-wave surveys are now becoming avail- 
able (Barger et al. 1998; Eales et al. 1998; Holland et al. 
1998; Hughes et al. 1998). 

The surface density of galaxies with flux densities 
greater than 4mJy at a wavelength of 850 /im is (2.5 ± 
1.4) x 10 3 deg -2 (Smail et al. 1997). Counts of lensed galax- 
ies that are consistent with these observations, which do 
not include the effects of lensing by discs, were presented 
by Blain (1998a,b,c). The local population of IRA S galaxies 
(Saunders et al. 1990) appears to undergo pure luminosity 
evolution of the form (1 + z) 7 where 7 ~ 4 out to z ~ 2.5 
(Blain et al. 1999a). This model is consistent with counts 



at other wavelengths in the millimetre/submillimetre wave- 
band (Wilner & Wright 1997; Kawara et al. 1998; Lagache 
et al. 1998), and with the intensity of diffuse background ra- 
diation in the submillimetre/far- infrared waveband (Puget 
et al. 1996; Guiderdoni et al. 1997; Schlegel, Finkbeiner & 
Davis 1998; Hauser et al. 1998; Fixsen et al. 1998). 

4.2 Calculating counts of lensed galaxies 

We determine counts of lensed galaxies by using equation (3) 
with $' calculated using F d i sc (equation 15), normalised as 
described above, in both models of disc evolution described 
in Table 2. The population of distant dusty galaxies assumed 
is described by the 'Gaussian' model of Blain et al. (1999a), 
which accounts for all the available submillimetre-wave and 
far-infrared data currently available. 

In order to illustrate the expected form of the counts of 
both lensed and unlensed galaxies at wavelengths of 1300, 
850, 450 and 60 /im, the counts of galaxies expected in the 
absence of discs are compared in Fig. 5(a). Lensed counts 
are calculated, assuming that M* does and does not evolve. 
If M* does not evolve, then the predicted count of lensed 
galaxies is greater. As a result, all previous predictions have 
been made assuming that M* evolves to produce conserva- 
tive estimates of the counts of lenses. To demonstrate the 
effects of including lensing by discs, the relative size of the 
counts of lensed galaxies with discs as compared with the 
counts shown in Fig. 5(a) are shown in Fig. 5(b). The addi- 
tional bias arising from disc lensing is typically a factor of 2 
and 3 in models A and B respectively. 

The only data currently available on the population of 
lensed dusty galaxies and quasars is provided by the identifi- 
cation of two luminous distant 60- /im IRAS sources as lenses 
(Close et al. 1995; Irwin et al. 1998; Lewis et al. 1998), and 
the identification of one known lensed quasar with a 60-pim 
IRAS source (Kneib et al. 1998). This indicates a 60-pim 
count of lenses of at least 7 x 10~ 5 deg~ 2 at a flux density of 
0.2 Jy, fully consistent with the value of about 7 x 10~ 4 deg -2 
shown in Fig. 5(a). 

The maximum magnification assumed in these calcula- 
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tions for a source at z = 1 is /i ma x(l) = 40, as assumed 
by Blain (1996). A larger value of /x max will lead to a more 
significant effect of lensing by discs, as discs are progres- 
sively more efficient at producing large magnifications as 
compared with SIS haloes: see Figs 1(a) and 4. This effect is 
emphasized in Fig. 6, in which the additional magnification 
bias arising from lensing by discs is compared for different 
values of /i max (l): 10, 20, 40 and 100. Observations made 
using millimetre-wave interferometers (Solomon et al. 1997; 
Mirabel et al. 1998; Sakamoto et al. 1999) indicate that the 
most luminous far-infrared emitting regions of distant galax- 
ies are rather small, with scales of less than several hundred 
parsecs. As /i max (l) — 25 is expected for a 1-kpc source 
(Peacock 1982), (1) = 40 seems to be very reasonable. 
Larger magnifications may also be possible, at least for any 
very compact subsample of distant submillimetre-luminous 
galaxies and quasars, such as the 'extreme starbursts' im- 
aged by Downes & Solomon (1998). In this case, the numbers 
of lenses with flux densities of order 100 mjy, and thus the 
magnification biases due to discs (Fig. 5b), will be increased. 

The fraction of lensed galaxies detected in a carefully 
designed flux-limited submillimetre-wave survey is expected 
to be several per cent (Blain 1996, 1998b), and it is reason- 
able to expect lensing by discs to increase this fraction by a 
factor of between 2 and 3. 



5 THE PROSPECTS FOR FUTURE 
OBSERVATIONS 

5.1 Potential instruments 

Many forthcoming instruments will be useful for carry- 
ing out submillimetre-wave lens surveys. Foremost amongst 
them are large ground-based millimetre/submillimetre inter- 
ferometer arrays (MIAs), such as the MMA (Brown 1996): 
see also Downes (1996) and Ishiguro et al. (1992). These 
instruments, with subarcsecond angular resolution, will be 
required in order to confirm submillimetre-selected lens 
candidates (Blain 1998b). Large single-antenna telescopes 
such as the US-Mexican 50-m Large Millimeter Telescope 
(LTM/GTM; Schloerb 1997) and a 10-m telescope at the 
South Pole (Stark et al. 1998), are also being designed and 
built. These telescopes will provide large and uniform focal 
planes, able to exploit large format bolometer array receivers 
(Glenn et al. 1998), and would be very useful for wide-field 
millimetre/submillimetre- wave surveys (Blain 1998c). The 
Planck Surveyor (Bersanelli et al. 1996) and Far-Infrared 
and Submillimetre Telescope (FIRST; Pilbratt 1997) space 
missions also have great potential. 

5.2 The detection rate of lenses 

The relative speed of detecting lenses using the instruments 
discussed above is considered in more detail elsewhere (Blain 
1998c). There is freedom to choose a survey strategy for all 
but Planck Surveyor, and, if optimized, detection rates of or- 
der 0.03-0.3 h -1 are expected. If the effect of lensing by discs 
is included, then these detection rates would be expected to 
increase by a factor of about 2 or 3. If the optical depth 
to submillimetre-wave lensing by discs is as large as we ex- 
pect, then follow-up observations of detected lenses should 




Limiting Flux Density / Jy 

(b) 



Solid lines: Model A Dashed lines: Model B 

Wavelength 60. 450. 850 and 1300/xm ^ ^ ^ 




10 4 10 3 0.01 0.1 1 

Limiting Flux Density / Jy 

Figure 5. (a) Counts of lensed galaxies expected in the millime- 
tre/submillimetre and far-infrared wavebands (Blain et al. 1999a). 
No effects of lensing by discs are included. Lensed counts with 
and without an evolving population of haloes, as parametrized 
by M*, are shown, (b) The additional magnification bias factor 
that should multiply the counts in (a) when lensing by discs is 
considered in both models of disc evolution (Table 2 and Fig. 4). 
Some curves stop before the right-hand edge of the figures because 
the predicted counts fall below 2.4 X 10~ 5 dcg~ 2 (l/47rsr _1 ). 

reveal that a large fraction are lensed by edge-on discs. 
Submillimetre-wave point sources will be extracted from the 
all-sky Planck Surveyor survey (Blain 1998b), down to the 
confusion limit (Blain et al. 1998a,b). More than a thousand 
lensed galaxies and quasars were predicted, without consid- 
ering the effects of lensing discs. This number should also 
be increased by a factor of 2 or 3 when lensing by discs is 
taken into account. 

5.3 Imaging lens candidates using the MMA 

Multiple lensed images formed by galaxies are separated on 
arcsecond scales. A large interferometer array, such as the 
MMA, would thus be the ideal telescope with which to inves- 
tigate lensing by distant galaxies in detail. In Fig. 7 we show 
simulated images of a 10-arcsec wide field, contained within 
the primary beam of the MMA, at an angular resolution of 
0.1 arcsec. The contour levels correspond to three times the 
850- /im la sensitivity limits in MMA integrations lasting 
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Figure 7. Simulated MMA images obtained at a wavelength of 850 um for a representative region of sky that is chosen to contain 
a strongly lensed source. About 2 per cent of fields centred on a foreground galaxy would be expected to contain such a background 
source. The lensed images produced by a Milky Way galaxy at z = 0.3, inclined 5° from the edge-on aspect, and the emission 
from the disc are shown in (a). In (b) the corresponding images produced by an SIS lens with the same mass as the disc, also at 
z = 0.3 are compared. The background population of sources are shown in (c). The contour levels correspond to the la noise levels 
expected in 1-min, 1-h and 24-h integrations: 220, 28 and 5.8 /iJy beam -1 respectively. 
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Figure 6. Differences in the additional magnification bias due to 
lensing by discs, expected at a wavelength of 850 (im in disc evo- 
lution Model B: see Table 2. A larger value of /i max corresponds to 
a more significant effect, emphasizing the results shown in Fig. 4. 

lmin, lh and Id: 220, 28 and 5.8 /iJy beam" 1 respectively. 
Note that the primary beam area of the MMA at 850 fim is 
about 3x 10~ 5 deg -2 , and so the brightest source in this area 
would be expected to have a flux density of 40 /iJy (Fig. 5a) . 
In Fig. 7(c) a typical region of the background sky is shown. 
The source distribution in this case is chosen so that one 
source will be magnified significantly by a disc galaxy at the 
centre of the figure. All sources are Gaussian and 0.2 arcsec 
across at full width at half- maximum. An alignment between 
source and lens this close or closer is expected in about 2 per 
cent of cases. In Fig. 7(a) the effects of lensing these sources 
by an edge-on disc galaxy, similar to the Milky Way, at a 
redshift of 0.3 is shown. The submillimetre-wave emission 
from the lensing galaxy is calculated by transforming the 
observation of a nearby edge-on disc galaxy (Israel, van der 
Werf & Tilanus 1998) to z = 0.3, giving a total flux density 
of about 1 mjy and a scalelength of 1 arcsec. In Fig. 7(b) the 



much less significant effects of lensing by an SIS dark matter 
halo with the same mass at the same redshift are compared. 

In addition to following up the sources detected in sur- 
veys using single-antenna telescopes (Blain 1998c), a deep 
targeted MIA survey of known distant edge-on discs, pro- 
vided for example by the Sloan Digital Sky Survey catalogue 
(Margon 1998), would be a very powerful way of detect- 
ing lensed images of very faint distant submillimetre-wave 
galaxies. Integration times of about 24 h per source would 
be required however. 



6 CONCLUSIONS 

(i) We have estimated the enhanced magnification bias 
expected for statistical lensing due to a population of bary- 
onic discs in dark matter haloes. The bias is expected to be 
enhanced by a factor of 2 or 3 as compared with a popu- 
lation of SIS lenses. This bias could be much larger if the 
sources are smaller than several hundred parsecs in size, and 
can thus be magnified by a factor in excess of 40. 

(ii) In the submillimetre waveband, lensed images will 
not suffer from the effects of dust absorption in the lens- 
ing galaxy. Hence images lensed by edge-on discs should be 
selected in an unbiased manner in submillimetre-wave sur- 
veys made using the MMA, LMT/GTM, a 10-m telescope 
at the South Pole, and the space-borne instruments Planck 
Surveyor and FIRST. 

(iii) The subarcsecond angular resolution of the MMA 
will be required in order to study lensed images in detail. 
A targeted search of known edge-on discs could provide a 
fruitful technique for detecting and observing any very faint, 
distant dusty galaxies. 
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